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Abstract 

We recently proposed a spontaneous A4 flavor symmetry breaking scheme imple- 
mented in a warped extra dimensional setup to explain the observed pattern of quark 
and lepton masses and mixings. The main features of this choice are the explanation 
of fermion mass hierarchies by wave function overlaps, the emergence of tribimaximal 
(TBM) neutrino mixing and zero quark mixing at the leading order and the absence 
of tree-level gauge mediated flavor violation. Quark mixing and deviations from TBM 
neutrino mixing are induced by the presence of bulk A4 flavons, which allow for "cross- 
brane" interactions and a "cross-talk" between the quark and neutrino sectors. 

In this work, we study the constraints associated with the recent measurements of 
# 13 9° by RENO and Daya Bay, forcing every model that predicts TBM neutrino 
mixing to account for the significant deviation of #13 from 0, while keeping the values 
of 612 and 023 close to their central experimental values. We then proceed to study in 
detail the RS-A4 contributions to \i — > e, 3e, generated at the tree level by virtue of 
anomalous Z couplings. These couplings arise from gauge and fermionic KK mixing 
effects after electroweak symmetry breaking. Since the experimental sensitivity for 
BR(ji — > e, 3e) is expected to increase by five orders of magnitude within the next 
decade, it is shown that the RS-A4 lepton sector can be significantly constrained. 
Finally, we show that when "cross-brane" interactions are turned off, the Z couplings 
are protected against all anomalous contributions and a strong correlation between #13 
and the deviation from maximality of O23 is found. 



1 Introduction 



Recently we have proposed a model |T] based on a bulk A A flavor symmetry [2j[3] in warped 
geometry [4], in an attempt to account for the hierarchical charged fermion masses, the 
hierarchical mixing pattern in the quark sector and the large mixing angles and mild hierarchy 
of masses in the neutrino sector. In analogy with a previous RS realizations of A4 for the 
lepton sector [5j[6] , the three generations of left-handed (LH) quark doublets are unified into 
a triplet of A4, this assignment forbids tree level FCNCs driven by the exchange of KK 
gauge bosons, as long as only brane localized effects are considered. The scalar sector of the 
RS-A4 model consists of two bulk flavon fields, in addition to a bulk Higgs field. The bulk 
flavons transform as triplets of A4, and allow for a complete "cross-talk" [7] between the 
A4 — > Z 2 spontaneous symmetry breaking (SSB) pattern associated with the heavy neutrino 
sector - with scalar mediator peaked towards the UV brane - and the A 4 — > Z 3 SSB pattern 
associated with the quark and charged lepton sectors - with scalar mediator peaked towards 
the IR brane - allowing us to obtain realistic mixing angles in the quark sector and deviations 
from TBM mixing in the neutrino sector. A bulk custodial symmetry, broken differently at 
the two branes [8] , guarantees the suppression of large contributions to electroweak precision 
observables p], such as the Peskin-Takeuchi S, T parameters. However, the mixing between 
zero modes of the 5D theory and their Kaluza-Klein (KK) excitations - after 4D reduction - 
may still cause significant new physics (NP) contributions to SM suppressed flavor changing 
neutral current (FCNC) processes. 



In 10 we have performed a thorough study of the RS-A4 contributions to one loop 
AF = 1 dipole transitions (nEDM, e'/ex, b — > 5(^)7) and tree level Higgs mediated FCNC 
AF = 2 transitions (K° — K° and D° — D°mixing). The main difference between the 
RS-A4 setup and an anarchic RS flavor scheme [II] lies in the degeneracy of fermionic LH 
bulk mass parameters, which implies the universality of LH zero mode profiles and hence 
forbids gauge mediated FCNC processes at tree level, including the KK gluon exchange 
contribution to ex- The latter provides the most stringent constraint on flavor anarchic 
models, together with the neutron EDM (nEDM) [TT ,12 . However, the choice of the common 
LH bulk mass parameter, is strongly constrained by the matching of the top quark mass 
(mt(1.8TeV)~ 140 TeV) and the perturbativity bound of the 5D top Yukawa coupling, y t . 
Most importantly, when considering the tree level corrections to the Zbb coupling against 
the stringent electroweak precision measurements (EWPM) at the Z pole, we showed [I] 
that for an IR scale, Am ~ 1.8 TeV and w 125 GeV, is constrained to be larger than 



0.34. Assigning = 0.4 and matching with m t we obtain \y t \ < 3, which easily satisfies 



the 5D Yukawa perturbativity bound. For Am ~ 1.8 TeV the lightest fermionic KK mode is 
b^\ the "fake" (cusodial) SU(2)r partner of the 5D fermion, whose zero mode is identified 
with tn, with M KK ~ 0.95 TeV. The second most significant constraint on the KK mass 
scale comes from one-loop Higgs mediated dipole operator contributions to the b — > 37 
process, Mj^f ~ 2.55A^ S7 > 1.3 Y TeV, where Y is the overall scale of the dimensionless 
5D Yukawa coefficients. The constraints on M^k coming from e'/ex (nEDM) were shown 
to be weaker by at least factor of 2(6) (See [To]). 
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In this paper we focus on the charged lepton and neutrino sectors of the RS-A 4 setup. 
We start by performing a systematic study of the RS-A 4 predictions for the neutrino mixing 
angles in light of the recent measurements of #13 ~ 9° by the RENO 13 , Daya 
and DOUBLECHOOZ [i~5] experiments, supporting previous indications in T2K 
MINOS " 
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16 and 



It is shown that significant deviations from TBM neutrino mixing are generated 
by both "cross talk" interactions in the charged lepton sector and higher order corrections 
to the heavy Majorana and (light) Dirac mass matrices. The model predictions are then 

Subsequently, we inspect 



18,19 



compared with the results of the most recent global fits 
the RS-A 4 contributions to the charged lepton flavor violating processes /i — > 3e and fi — > e 
conversion. We show that the dominant contributions come from tree-level Z boson exchange 
diagrams, where the flavor violating Z boson couplings are induced by gauge and fermion 



KK mixing effects. We then compare our predictions with the existing upper bounds 20 21 



and those expected to come from a series of experiments, currently under construction in 
J-PARC |22}|M], Osaka (RNCP) [25], Fermilab and PSI 

The paper is organized as follows. In Sec. |2l 



within the next decade. 



we summarize the most relevant features 
of the RS-A 4 setup needed for the analysis of the following sections. In Sec. [3] we study 
the modifications of the PMNS matrix induced by higher order effects in both the neutrino 
and the charged lepton sectors and perform a numerical scan to test our predictions against 
the experimental bounds. In Sec. |3.1 we specialize to the brane localized RS-A 4 setup. In 
Sec. [4] we study the effect of gauge boson and KK fermion mixing on Z couplings and obtain 



analytical estimations of the Zfxe coupling relevant for the /1 — > e, 3e processes. In Sec. 4.4 



we obtain the exact results by a numerical diagonalization of the zero+KK mass matrices 
and perform a scan over the parameter space of the RS-A 4 setup including/excluding cross- 
brane effects. It is shown that the brane localized version of RS-A 4 stays protected from KK 
mixing contributions to off diagonal Z couplings. We conclude in Sec.[5j 



2 The RS-A 4 model 



The RS-A 4 setup [T] is illustrated in Fig. [TJ The bulk geometry is that of a slice of AdS$ 
compactified on an orbifold S1/Z2 |I] and is described by the metric on the bottom of 
Fig. [l] All 5D fermionic fields propagate in the bulk and transform under the following 
representations of (SU(3) C x SU(2) L x SU(2) R x U(1) B -l) xA 4 xZ 2 |: 

e R © e' R © e" R ~ (1, 1, 2, -1) (1 © V © 1") (+), t L ~ (1, 2, 1, -1) (3) (-), 

u R © u' R © u" R ~ (3, 1, 2, §) (1 © V © 1") (+), u R ~ (1, 1, 2, 0) (3) (-), (1) 

d R © d' R © d" R ~ (3, 1, 2, §) (1 © V © 1") (+), Q L ~ (3, 2, 1, |) (3) (-). 

The SM fermions, including right handed (RH) neutrinos, are identified with the zero modes 
of the 5D fermions above. The zero (and KK) mode profiles are determined by the bulk mass 
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y = 



SU(2) L x U(1) Y 




Si/Z: 



a 4 x (z 2 y 




SU(2) t x SU{2)r x U(l) B-L 



ds 2 = e 2k ^j] ilu dx^dy v + dy : 



y — ttR 



SXJ(2) D xU(l) B - L 



Figure 1: A pictorial description of the RS-A4 setup. The bulk geometry is described by the metric 
at the bottom and k ~ Mpi is the AdS^ curvature scale. All fields propagate in the bulk and the 
UV(IR) peaked nature of the heavy RH neutrinos, the Higgs field, the t quark and the A4 flavons, 
<J? and Xi is emphasized. The SSB patterns of the bulk symmetries on the UV and IR branes are 
specified on the side (for A4) and on the bottom (for SU{2)i x SU{2)r x U(1)b-l) of each brane. 



of the corresponding 5D fermion, denoted by Cq Lm ^ itej k and boundary conditions (BC) (l0| . 
The scalar sector contains the IR peaked Higgs field and the UV and IR peaked flavons, x 
and $, respectively. They transform as: 



1,M, 



0) (3) 



X 



;i,l,l,0)(3)(+), H (1,2, 2,0) (!)(+). 



(2) 



The SM Higgs field is identified with the first KK mode of H. All fermionic zero modes 
acquire masses through Yukawa interactions with the Higgs field and the A4 flavons after SSB. 
The 5D (SU(3) C x SU(2) L x SU(2) R x U{1) B -l) X A A x Z 2 invariant Yukawa Lagrangian 
will consist of only UV/IR peaked interactions at the leading order (LO), while at the next 
to leading order (NLO) it will also give rise to "cross-talk" and " cross-brane" effects [l]. 
The LO interactions in the neutrino sector are shown in [l] using the see-saw I mechanism, 
to induce a tribimaximal (TBM) [29] pattern for neutrino mixing. Below we are going to 
show that "cross brane" and "cross talk" interactions, induce deviations from TBM, such 
that the model predictions for #12,23,13 are in good agreement with the current experimental 
bounds 18,19 . Starting from the quark and charged lepton sectors, the most relevant terms 



of the 5D Yukawa Lagrangian are of the following form: 
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5D D k 2 \Ql^l)^U[u r ,d R ,e R )H (Q L J L )®xti{u R ,a R ,e R , 



£lo Cnlo 

(3) 

Notice that the LO interactions are peaked towards the IR brane while the NLO interactions 
mediate between the two branes due to the presence of both $ and x- 
The VEV and physical profiles for the bulk scalars are obtained by solving the corresponding 
equations of motion with a UV/IR localized quartic potential term and an IR/UV localized 
mass term 30 . In this way one can obtain either UV or IR peaked and also flat profiles 
depending on the bulk mass and the choice of boundary conditions. The resulting VEV 
profiles of the RS-A 4 scalar sector are: 

v %> 9) = F o (0o)e (2+fe W )fe(l2/| - 7rii) vf = Xoe {2 -^ )k]v] (l - e Wx)H\v\-*R)) , (4) 



where 0h,<s>, x = y 4 + $ x and Hh,$,x is the bulk mass of the corresponding scalar in units 
of k, the cutoff of the 5D theory. The following vacua for the Higgs and the A± flavons $ 
and x, 

($) = {v h v h V4> ) (x) = (0, v x , 0) (H) = v H ( \ \j , (5) 

provide at LO zero quark mixing and TBM neutrino mixing (2j[7j. The stability of the above 
vacuum alignment is discussed in [I]. The VEV of $ induces an A 4 — » Z 3 SSB pattern, which 
in turn induces no quark mixing and is peaked towards the IR brane. Similarly, the VEV 
of x induces an A 4 — > Z 2 SSB pattern peaked towards the UV brane. Subsequently, NLO 
interactions, including both $ and x-, break A 4 completely and induce quark mixing and 
deviations from TBM. The Higgs VEV is in charge of the SSB pattern SU(2) L x SU(2) R -> 
SU(2) D , which is peaked towards the IR brane. The (gauge) SSB pattern on the UV brane 
is driven by orbifold BC and a Planckian UV localized VEV, which is effectively decoupled 
from the model (8). 

To summarize the implications of the NLO interactions in the quark and charged lepton 
sectors, we provide the structure of the LO+NLO mass matrices in the ZMA (l): 

( q AD AD AD \ / fu,d,eu,d/ fc,s,fj.u,d,t f t,b,T u,d,i 

Du,d,e i>c,s,ii Utfi.T \ I Jx 1 ■> X 2 J x 3 

vt% ^vlt + 

vi% <A c 4 ,t uy% ) \ f^ y ^ fwy%u f x ^yuA,t 

V v y 

V3U(u)diag(y* D d . e .) 

(6) 

where u = e 2m/3 , v = 174 GeV is the 4D Higgs VEV, yt% AsAe ^ r are the effective 4D 
LO Yukawa couplings, (xf,yf) = (x\,y^y A ® Ai Jy UiAi ^ and y Ui , di>ei are the 5D 

LO (NLO) Yukawa couplings. The function fgA* ~ 2/3 x C x /(12 - c Ql/l - c Ui4i , ei ) ~ 0.05 
accounts for the suppression of the NLO Yukawa interactions and C x = Xo/Mp 2 ~ 0.155. 



4 



Finally, the unitary matrix, U{uS) is the LO left diagonalization matrix in both the up 



and down sectors. 



)lOi which is independent of the LO Yukawa couplings, while 
iy^' a,l )LO — l(see |lj). Using standard perturbative techniques on the matrix in Eq. (|6 



we obtained (YlIr^nlo jl.10 at o(j. 
matrix is given by 



. The left-handed diagonalization 



f 



f^(xf + cuyf 



(7) 

where q = u,d. 

Similarly, the right diagonalization matrices in the quark and charged lepton sectors, to 
first order in f^ i,di,ei {x^ ,d,e ,y^' d,e ) are given by 



V, 



q ,e 



R 



1 



where q = u,d and the A^' are given by: 



A q ' e 



A 



q.( 



-(Af)* -{Aff 



A\ 



q,£ _ m u,d,< 



C,S,fl 



f u, d ,e (( x fy + co 2 (yfy)+f : 



A 9 '' 

^2 



A3,' 
^3 



■m u ,d,< 



u,d,e 
X 



~q,t\* 



+ u{y\ 



1 if + yt 



fi^^t + yl 



m 



c,s,n 



m,b,r 



The suppression by quark mass ratios of the off-diagonal elements in V^ ,d , stemming from 
the degeneracy of LH bulk masses, turned out to play an important role in relaxing the flavor 



(9) 

(10) 
(11) 



violation bounds on the KK mass scale, as compared to flavor anarchic frameworks 10 



3 The neutrino sector — Higher order corrections to the 
PMNS matrix and 13 

The global fits based on the recent measurements of — > v e appearance in the RENO, 
Daya Bay, T2K, MINOS and other experiments, allow one to obtain a significance of 10<r 
for 6*i3 > 0, with best fit points at around #13 ~ 0.15, depending on the precise treatment of 



reactor fluxes 18 , 19 . We wish the RS-A4 higher order corrections to the PMNS matrix to 
be such that the new fits are still "accessible" by a significant portion of the model parameter 
space. 
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In the RS-A4 model deviations from TBM neutrino mixing are generated by higher order 
effects on the UV/IR branes and by "cross-talk" and "cross-brane" effects. In [l] we studied 
in detail the textures induced by the three types of higher order corrections. Here we 
quote the dominant corrections to the Dirac and Majorana mass matrices and the resulting 
modification of V£ the left diagonalization matrix of the effective (LL) neutrino Majorana 
mass matrix. 

Starting from the Dirac mass terms we have: 



C D V C k- l ' 2 yJ L Hv R + k-WhHxvR + k- 7 / 2 y% 2 e L H$ 2 v R , 



'121 



where y u , y^ and y* 2 are dimensionless 5D Yukawa couplings. The 4D effective mass matrix is 
obtained by integrating over y the profiles of the fields/VEV's in each of the above operators 
and acquires the following form [I]: 



D 



m 



D 



1 + ei e 2 e 3 + e x 
e 3 1 + ei e 2 



(13) 



where the m^, e x and 61,2,3 entries come from the first, second and third terms of Eq. (12), 
respectively. Notice that the smallness of e x , 1,2,3 compared to comes not only from the 
(mild) suppression of the flavon VEV's with respect to the 5D scale k ~ O(Mpi) but also 
from the associated overlap correction factors [1,10 . 



Turning to the Majorana mass terms we have: 

C A J C Mv%v R + k- x l 2 y x xv R v R + k- 2 y x 2 X 2 u R u R , (14) 
The resulting 4D (heavy) Majorana mass matrix acquires the following form: 



M, 



M 




(15) 



where M, M x and €4,5 come from the first, second and third terms of Eq. (14). The smallness 



of €4,5 comes mainly from the suppression associated with the flavon VEV's, since the effect 
of overlap corrections is milder for UV peaked operators fil. 

The effective (light) Majorana mass matrix is now obtained using the see-saw mechanism, 
Mfff = -(M®) T (MM)- l M®. If we turn off ei_ 5 , x we obtain the LO results of HI, for which 
the mass spectrum is given by 



ffjdiag. 



{v v L ) T M: ff vi 



D\2 



M 



-diag 



l + q' 



(16) 



where q = M x /M. Notice that the overall neutrino mass scale is set by the ratio (m®) 2 /M, 
while the type of hierarchy and its extent are determined by q. Using Am 2 tm ~ 2.43- 10 _3 eV 2 



6 



and Am^ oi ~ 



7.6 • 10~ 5 eV 2 from [18] , we are able to obtain a cubic equation for q with 
solutions q ~ [—2.02,-1.98,0.79,1.19], where the first (last) two solutions correspond to 
inverted (normal) mass hierarchy llj. The LO diagonalization matrix V£ is given by: 



vr 



l/y/2 -l/v 7 ^ 
1 

l/y/2 l/y/2 



(17) 



We now proceed to obtain the corrections to V£ considering the higher order effects 

0.05, we can perform (analytically) 



encoded in Eqs. (13) and (15). Since numerically e% 



-5,x 



the perturbative diagonalization of {Ml 



eff. 



NLO 



to leading order in each of the e's. The 



resulting (Vj/ 



NLO 



VI is given by: 



/ J_ _l e 2- £ ; 
' V2 + 2V2 



\V2 



'2V2q 



2 -"(« • < 3 

'2V2q 



2V2 



+ i 



Im(e 4 ) 



Im(e 4 ) 



£2+^3 
9 



£2+<;3 
1 



+ ilm(e 2 ) + Re(e 3 ) 
1 

+ ilm(e 3 ) + Re(e 2 ) 



1 +^r+^Im(e 4 )\ 



V2 



2V2 2V2q 



Using V£ of Eq. (|7j) and (V^) of Eq. (18), we obtain the modified PMNS matrix, which 
includes all dominant higher order effects from the charged lepton and neutrino sectors, 
Vpmns = W = (VffiV£. Notice that all elements of the PMNS matrix will be still in- 
dependent of the LO Yukawa couplings (y e ,n,T,v,x,M), due to the Form diagonalizability of 
the LO mass matrices. We avoid writing the resulting PMNS matrix explicitly, due to its 
complexity. By expanding the basis independent relations 

13 = arcsm|W 13 |, 6 12 = axctan(|Wi 2 |/|Wii|), 9 23 = a icsm(\W 23 \/\W 33 \), (19) 

we are able to obtain approximate analytical expressions for $12,13,23, 



sin 9 



12 



1 + 



2 + q 
3q 



e 2 + e 3 + 2Re(e 2 + e 3 ) ) - ( co 2 (x e 2 + y e 2 ) + u{x\ + y\ 



(20) 



sin 9 



13 



y/6 



J Y 



(4 + yi)(l-u,) + (4 + yi)(l-u, 2 ) 



, ^ 2 ~1) T 
H 111164 

Q 



(21) 



sin 9- 



2.3 



1 

71 



—uj 



e 2 - e 3 J 



J x 



uj 2 x 3 + wzf + 2Re(yl\ 



i(q 2 - 1) 



Im e 4 



• (22) 
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Figure 2: Model predictions for 823 vs. 813 for normal (left) and inverted (right) mass 
hierarchy, including all dominant higher order and cross talk effects. The purple points 
represent the entire sample, while the red points satisfy the 3a bounds on 8\2- The rectangles, 
confined between the horizontal and vertical dashed green lines, represent the 3a allowed 
regions from the global fits of i 18, 19l 



We are now ready to perform a systematic study of the RS-A 4 predictions for #12,13,23- 
To this end we recall that ei-5, x of Eqs. (13) and (15), implicitly contain the dimensionless 
Yukawa couplings of their corresponding operator in Eqs. ( 12 ) and ( 14 ) . For this reason ei_5 iX 
will be multiplied by an 0(1) complex Yukawa coupling for the purpose of the numerical scan 
below. We generate sample of 20000 points, in which all Yukawas (xi,yf,y ei ) are complex 
numbers with random phases and magnitudes normally distributed around 1 with standard 

Using the overlap integrals performed in (T, 10 



deviation, a = 0.5. 



we get e x 
in Eq. 



0.07, 

we 



ei,2,3 ^ 0.04, e 4 , 5 ^ 0.09 and f^ T = V3f^' T ~ 0.08. By Substituting ei_ 5 , x 
study the correlations between the various neutrino mixing angles. The results of tne scan 
are plotted in Figs. [2]-[4| The rectangle confined between the dashed green lines represent the 
3a allowed regions from the global fits of 18 19 . The purple points in each plot represent the 



entire sample, while the red points, depicted on top of the purple ones, satisfy in addition the 
3a constraint from the mixing angle not appearing in the same figure. In other words, the red 
points in Fig. [2] satisfy the 3a constraint on # 12 and so on. The best fit point (s) are depicted 
in black (blue) for (#23)best values from the first (second) octant, where A^ack < ^Xbiue i n 
the global fit of (18). 

We realize that there are only mild differences between the NH and IH cases, coming 
from the terms proportional to (2 + q) and (q 2 — 1) in Eqs. (20) -(22). Qualitatively, this 
means that the deviations from TBM mixing in the IH case are slightly larger for #13,23 and 
slightly smaller for # 12 , as can be seen in Figs. [2]-[4| To quantify the viability of the RS-A 4 



setup in predicting realistic neutrino mixing angles, we can define a success rate for the NH 
and IH cases in analogy with 31 . Namely, we ask ourselves what is the percentage of the 
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Figure 3: Model predictions for # 12 vs. 813 for normal (left) and inverted (right) mass 
hierarchy, including all dominant higher order and cross talk effects. The purple points 
represent the entire sample, while the red points satisfy the 3cr bounds on 823- The rectangles, 
confined between the horizontal and vertical dashed green lines, represent the 3a allowed 



regions from the global fits of 1 18, 19l 



points, £, satisfying the 3a constraints for all neutrino mixing angles. We get £nh — 10% 



and £ih — 12.5%, which is slightly better than the results for typical A 4 models in 31 . We 



avoid analyzing the RS-A 4 predictions for 5q P , due to the lack of sufficient experimental 
data and its potential sensitivity to higher order effects of 0(ef,f^). The above analysis re- 
demonstrates |32}[33| the viability of models predicting TBM at LO, with corrections coming 
from higher order effects in both the charged lepton and neutrino sectors and in particular 
the RS-A 4 setup. 



3.1 Simplifications in the brane localized RS-A4 setup 



If we specialize to the brane localized version of the RS-A 4 model |5j|6] the cross talk inter- 
actions associated with the x l 2 3 and y^ 3 are absent and the expressions for the mixing angles 
simplify to the following form: 



sin^ 



13 



itf ~ 1) 

V6q 



sin 8 



23 



Im e 4 



2^2 



sin 8 



12 



+ i 



2V2 



+ 



3V3q V 



£2 - 63) + 




(23) 
(24) 
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Figure 4: Model predictions for # 2 3 vs. Qyi f or normal (left) and inverted (right) mass 
hierarchy, including all dominant higher order and cross talk effects. The purple points 
represent the entire sample, while the red points satisfy the 3cr bounds on #i 3 . The rectangles, 
confined between the horizontal and vertical dashed green lines, represent the 3a allowed 



regions from the global fits of 1 18, 19l 



We immediately realize that #13 is now governed by q and Im(e 4 ). We take e 4 to be pure 
imaginary and fix its magnitude to obtain the best fit value sin 2 (#i3) ~ 0.158 18 19 . The 



fixing of e 4 depends on the type of hierarchy, as determined by q = M x /M (Eq. (16)). For 
the NH case we get Im(e4) ~ 0.5(0.97) for q ~ 0.79 (1.19), while for the IH case we have 
Im(e 4 ) ~ 0.6 for q = —2.02 (—1.98). These solutions for Im(e 4 ) implicitly include y ei , the 
dimensionless 5D Yukawa coupling of the operator k~ 2 x 2 ^R^R, entering as e 4 = k~ x y^v^ in 
the 4D effective theory on the UV brane. For the characteristic value k~ 1 v^ v ~ 0.1 we have 
\y^ H \ — 5(10.5) and \y^\ — 6, all of which satisfy the naive dimensional analysis (NDA) 
perturbativity bounds from [5] |y £4 | < Air. 

Now that we solved for 13 , the values of #12,23 are controlled by the e 2 ,3 parameters 
coming from the operator k~ 2 ii,H<& 2 VR. In general, we want #23 to deviate towards its best 



fit value in the first octant sin 2 9^ u ~ 0.427 19 and #12 stay close to its trimaximal value. 



Observing Eq. (24), we realize that #13 and #23 contain exactly the same correction term 
proportional to Im(e 4 ), fixed by the value of #13. Thus, we have a strong correlation between 
the increase of #13 from zero and the deviation of #23 towards the first (or second) octant. 
Since £2,3 are of the same strength, we can take them to be pure imaginary and opposite in 
sign to change #23 to the desired value (623 s * 7 ) while keeping sin 12 — l/v3. Since the latter 
possibility is highly fine tuned, we find it more instructive to consider the more general 
case in which 62,3 — 0.04 y €23 , where y €23 are complex numbers with random phases and 
magnitudes normally distributed around 1 with standard deviation of 0.5. In this case, #12,23 
are governed by four parameters, two phases and two Yukawas. We first specify to the case 
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Figure 5: The brane localized RS-A4 predictions for # 12 vs. 623 for normal (left) and in- 
verted (right) mass hierarchy, including all dominant higher order effects. The best fit value 
sin(#i3) ~ 0.158 fixes Im(e4) up to a sign. The purple (red) points correspond to a positive 
(negative) value of Im(e4) and we scan over the phases of €2,3 fixing their magnitudes to 
0.04 to make the comparison with the "cross-talk" case more transparent. The white rectan- 
gles, confined between the horizontal and vertical dashed green lines, represent the 3a allowed 



regions from the global fits of 1 18, 19l 



in which the magnitudes are fixed, where the number of observables equals the number of 
input parameters, in order to look for a pattern (see also [34]). 

We depict the results in Fig. [5j It can be seen that 6*23 moves strongly to the first (second) 
octant if Im(e4) is positive (negative) and that the best fit point is rather far from the center 
of the resulting distributions. It is not sensible to talk about a success rate in such a case 
due to the assignment imposed by #13 and the simplifying assumptions on £2,3. Relaxing 
the assumption on the magnitudes of £2,3 have the effect of spoiling the tendency towards 
the first or second octant, as can be seen by the slight overlap between the "positive" and 
"negative" cases, as can be seen in Fig. [6j 

The analysis of this section shows that the brane localized version of RS-A4 is more 
constrained by the experimental neutrino data, due to the fact that the Yukawa coupling 
associated with e 4 is fixed by 9 13 to values of order \y e4 \ ps 6 which are still below the 
perturbativity bound implied by NDA but still render the 5D theory itself less accurate. On 
the other hand, the small number of parameters (62,3,4) implies significant simplifications and 
an interesting correlation between the value of 6*13 and the (bi-directional) deviations from 
maximality of 623. This feature is attributed to the A4 assignments and vacuum alignment. 
Finally, the corrections to 612 in the IH case are extremely small due to the absence of 64 
and the (2 + q)/q term in Eq. (20). 
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Figure 6: The brane localized RS-A 4 predictions for 8 12 vs. 6 2 3 for normal (left) and inverted 
(right) mass hierarchy, scanning over the phases and magnitudes of €2,3. The best fit value 
sin(#i3) ~ 0.158 fixes 1111(64) up to a sign. The purple (red) points correspond to positive 
(negative) values oflm.(e±). The white rectangles, confined between the horizontal and vertical 
dashed green lines, represent the 3cr allowed regions from the global fit of [. 




4 The charged lepton sector- Anomalous Z couplings 
and cLFV 

Charged lepton flavor violating processes like /x — > 3e and // — > e conversion in the presence 
of a nucleus are the main focus of a series of experiments, planned to be performed in J- 
Osaka (RNCP) (25] 



PARC 22 24 



26,27 and PSI 28 , within the next decade 



Fermilab 

and increase the accuracy for these searches by roughly five orders of magnitude. The current 
upper bounds were set by the SINDRUM (ji -»■ 3e) |20) and SINDRUM II (// -»> e in Tig) 
21 experiments. 



As stated in the introduction the main source of charged lepton flavor violation (cLFV) 
in RS-A4 are anomalous Z couplings, induced by gauge and fermionic KK mixing, generating 
Tree level Z exchange contributions to \i — > e, 3e, B s — > ji + n~ and other processes. The effect 
of EWSB on the mixing of the Z boson with its KK partners and those of the (custodial) Z', 
can be studied directly in the zero mode approximation (ZMA), by solving the Z equations 
of motion in the vicinity of an IR boundary term which induces the SSB pattern SU(2)l X 
SU(2)r — > SU(2)d (e.g. [35]). On the other hand, the effect of KK fermion mixing has to 



be studied by considering the full KK mass matrices 36 . 



4.1 The impact of gauge boson mixing on Z couplings 

The effect of EWSB on gauge bosons mixing translates via the EOM to a distortion of the 
Z wave function near the IR brane . The advantage in using this procedure is that all gauge 
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bosons and their KK partners are automatically given in the physical basis. The canonically 
normalized physical Z boson wave function (profile) is then given by [35] 



f (o) 
J z 



(y) 



'txR 



1 + 



Mle 2ky 



Ak 2 



1 - 2ky 



(25) 



The modified Z wave function will induce small deviations in the Z coupling matri- 
ces A L>R (Z) } due to the different overlaps of the zero and KK modes profiles with /j (y). 
Schematically the overlap integrals associated with the Z couplings are written as: 



{ah 



'ttR 



{9l,r)smRui 
(9l,r)sm 



dyff 



l x~ (0,1,1 -- ■••)/ \ -(0,1,1- 

(y)XL,R {Ce,(e^T),y)X L ,R 



where (gf tR )sM = (^ 3 )< W _( ?L sin2 @w and X^r( c ^,t, y) are the canonically normalized 
wave functions for the zero mode fermions and their KK partners (denoted by m and n). 
Before EWSB the Z profile is flat and as a consequence the overlap integrals i?f m reduce 



,!+-,!-+. ..) 



Q,(e,/^,r) 



y), 

(26) 



to an integral over the n th and m th KK modes, satisfying R^ m (fz = l/VftR) = 5 nm by 
virtue of the orthonormality of the wave functions and implying [^^(Z)] oc lQ When 



EWSB takes place, is given by Eq. (25) and the coupling matrices, [y4 ifi (Z)] are no 
longer proportional to the identity As a result, once the Z coupling matrices will be rotated 
to the physical basis of the zero+KK fermions, new off diagonal entries will be generated in 
Al,r(Z), corresponding to flavor and KK number (and BC type) violating Z couplings. 

We first focus on the gauge boson contribution to flavor violating Z couplings in the 
ZMA, since in this case we are able to obtain simple analytical expressions using the left and 
right diagonalization matrices of Eqs. ^ and (|8]) - (JTTj). We start by performing the integral 
in Eq. (26) for the zero modes of £l and e^'"' ) (Eq. (Ill), which are identified with the charged 



lepton sector of the SM. Using 



(0) 'c,y) 



Xl,r 



fc(l -2c) 



,(l/2-c)%| 



-(0) / 

xl,r\ c i y 



■kR) = V2kF(c 



e fc7r_R(l-2c) _ I 

the resulting deviations in the diagonal entries of the Z coupling matrix are given by: 

(c) /, „ 5 -2c 
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zff 



^R(g z 



l,r)sm 



-"-00 — 



(9 



z \ 
l,r)sm 



M 2 Z F 2 ' 



A? fl (3 



knR- 



(27) 



(2f 



2c) V. 2(3- 2c) _ 

c e,e,fi,r- The physical couplings are obtained by rotating the Z 



where / = £, e, fi, r and c 

coupling matrices for the zero modes, [A L:R (Z)]^ Q to the mass eigenbasis: 



[^)] G o 



[Ar(Z)} 



(Vl)HA L (Z)} 



oo v l 

G yi 

00 V R 



^R(g z L )sM(V[y {R* (c e )-l)vi, 



00 ( 



(29) 
(30) 



1 At this stage we ignore the presence of opposite chirality ( ) KK fermions and (H — , — h) "custodians" , 

which will be thoroughly discussed in the next section. 



13 



Since [Al(Z)] 00 oc 1, no LH flavor violating Z couplings will be generated through gauge 
boson mixing after EWSB and in particular 8{g^) zauge = 0. On the other hand, the non 
degeneracy of RH bulk masses imply [v4 fi (Z)]^ ^cl and hence (Sg^) gaiUge 7^ 0. By using 
Eqs. (|8|-((9| and (28) -(30) we are able to estimate the strength of this effect 



« = (9 z / e : 



RS-Ai - {9 R , 



SM 



G 



J 12,21 



-5g 



R 



(31) 



where 5g R 1 ' J '' Zee are defined in an analogous way and extracted directly from Eq. (28). Since 
F(c e )/F(c IJ ) ~ rrie/m^, the dominant contribution to {5g R ) gaU ge will come from Sg^ 1 ^. 



We realize that 
characteristic values, A IR 



R ) gauge 



is suppressed by both F (c^M^/ Aj R , and f^[m e jm^). For the 
1.5 TeV, a ~ 0.51 and « 0.63, we have 8(g^) gauge « 5xl0- n . 



4.2 The impact of KK fermion mixing on Z couplings 

To account for KK mixing effects, we write the LO KK mass matrix for the first generation 
in the charged lepton sector in analogy with 10,12 . This mass matrix includes the zero and 



first level KK modes and acquires the following form: 10 

( 1 ' ^ ( y e FeF e r 00 x 
y e F e r 1Q x 



M 



KK 



KK, 



f(0) 
/e(l) 

g(l~) 



/ 



V 









1 

y* u r 222 x 



y e F e r ix 

Ve r 22X 

1 





VvTlllX 



1 / 



/ 



V 



'R 
e(l— ) 
R 

e (1) 

e R 



(32) 



/ 



where we factorized a common KK mass scale M KK ~ 2.45A/ R , y e>v = 2y e ^v$ D / A IR and 
the perturbative expansion parameter is defined as x = v/M KK . In the above equation the 
various r's denote the ratio of the bulk and IR localized effective couplings of the modes 
corresponding to the matrix element in question. For simplicity, we define rm = r n -+, 
rioi = r oi-+> r 22_= r i-i-i r 222 = ?"i-i+- and the notation for the rest of the overlaps is 
straightforward [10] . The associated Yukawa matrix, Y KKl is obtained by simply eliminat- 
ing x and the l's from the above matrix and once rotated to the mass basis it describes 

the physical coupling between (++), ( ), (H — ) and ( — h) KK modes. In 10 we have 

performed an analytical diagonalization of all one-generation zero+KK mass matrices in the 
quark sector and showed that the mild differences between the generations, stem from the 
structure of the 5D Yukawa couplings and non degeneracy of the RH bulk mass parameters. 

The full three generation mass matrix will be 12 x 12 and of similar structure, which is 
modified mainly by the A 4 flavor structure. We first decompose the 12 x 12 mass matrix 
of the RS-A4 charged lepton sector in terms of the one-generation KK+zero matrices in 



Eq. (32) 



M? K /M KK 



M 



Full 



M, 



KK 



xY, 



{ti?,F P 



KK\U2\ 5 

xY* K (yi°,F e ) 



X ^KK {Vl2°i Fp.) 

M? K /M j 



KK 



X ^KK (2/32 5 F/j, 



x \kk{Vi^ Ft 

x YkK (.V23 5 F T 

M? K /M KK 



(33) 
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where the expression in parenthesis for each off-diagonal element denotes the replacements 
to be made in Eq. (32) and the Mf^f matrices are of exactly the same structure in Eq. ( 32b . 
To account for NlAJ Yukawa interactions, the replacement yfP — » yff + y^ LO for the 
dimensionless 5D Yukawa coupling matrices applies. In the above equation, M^k is the KK 
mass corresponding to the degenerate left-handed bulk mass parameter, cf, and we normalize 
all matrices accordingly while still keeping the small non- degeneracy of the KK masses, 
M KK = ~ ~ 0.87Mfi ~ 0.95M^ ~ 1.08Mj£. Finally, v% D /A IR ~ 0.27 

and we assign q ~ c r ~ 0.51, c e ~ 0.787, c M ~ 0.63, y e ~ 1.2, ~ 1.55 and y T ~ 1.7 to 
reproduce the charged lepton masses with an IR scale A IR ~ 1.5 TeV. The resulting masses 
are m e = 0.511 MeV, m M = 105.7MeV and m T = 1.77 GeV. 

The reason fermion KK mixing is so important for off diagonal Z couplings is the presence 

of opposite chirality ( ) KK states and "fake" custodial partners, which are the SU(2)r 

partners of vr (Eq. ([!])) with ( — h) and (H — ) boundary conditions. Recall that a 5D fermion 
corresponds to two chiral fermions in 4D. As a result we will have LH states which couple 
to the Z boson as RH states (e.g. £^ ') and vice versa (e.g. e^ + ~^). As a result, the Z 
coupling matrices, Al,r{Z), which are proportional to the identity matrix in the ZMA (up 
to small deviations induced by Eq. (28)), contain instead both types of (diagonal) entries 
{9l)sm, (9r)sm and will thus acquire non vanishing off diagonal elements, once rotated to 
the common mass basis of the KK and zero mode fermions. 

Unfortunately, the three generation charged lepton (KK+zero) mass matrix can only be 
diagonalized numerically, due to its dimension and the large number of input parameters. 
The numerical diagonalization of this matrix is quite involved due to the presence of large 
hierarchies, ranging from m e to Mkk- It is thus essential to obtain an analytical approx- 
imation for {$Q*Jjr)kK) to get a control over the numerical results. For this purpose, we 
follow [36] and use an effective field theory approach, where all heavy modes are integrated 
out. Within this approximation the zero mode mass matrix is modified in the following way, 



(M 



on ) 



\KK 



M 00 + M fc MjT 1 M fc j M~ 1 M, 



1 

2 L 



Mo k M~ 2 Ml k Moo + M 00 Ml M- 



kO 



+ 



(34) 

where the first correction originates from the pure heavy (KK) mass terms, while the sec- 
ond correction containing Moo is generated by the redefinitions of the light (SM) charged 
lepton fields, which are required to bring the kinetic terms into their canonical form 36 . 
Notice that M o is precisely the ZMA mass matrix given in Eq. (|6]), upon the identifica- 



tion y\ 



4D 



tye.UirFlFe, 



AD 



/A 



IR) r 00 • 



In Eq. (|34|) the zero+KK mass matrix of Eq. rt33h 

; is obtained by 
To simplify the 



is decomposed into its "generational" 3x3 building blocks, such that M^j is obtained 
taking the (kj) element of each of the 4x4 block matrices in Eq. (33) 



notation, we label the zero and KK modes in the basis vectors of Eq. (32) as (0,1,2,3) 
Observing Eq. (32), we immediately realize that M ± = M 2 o 
1 



M: 



Mi 



1,33 



KK 



30 



M 



23 



M 



32 



0. 



1) 



and M 22 



diag(M^ 



L KK 



The modified mass matrix of Eq.(34) induces corrections to the left and right handed Z 
coupling matrices for the zero modes, [Al,.r(^)]ooj where the first non vanishing corrections 
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can be thought of as coming from placing two KK-zero mass insertions on the fermionic legs 
of a Zfjj R fL R vertex. The corrected coupling matrices are given by 



36 



iMz)t K 



[MZ)\1 K 



[A L (Z)} 00 + M 0k M k 2 [A L (Z)} kk Ml 



Ok 



M 0fc M fc - 2 Mi [A L {Z)\ m + [A L (Z)} m M m M m M^Ml 



[A R (Z)] Q0 + Ml M- 2 [A R (Z)] kk M, 



kO 



Ml Q M- 2 M k0 [A R (Z)] 00 + [A R {Z)] m M m Ml M- 2 M } 



kO 



+ ... 



+ ... 



(35) 



(36) 



In the above equations, the corrected coupling matrices [Al,r(Z)]q Q k are given in the 
interaction basis and should be rotated to the physical basis by the usual transformation 
(Yl rY[^-l,r]oo K Vl r- To account solely for KK mixing effects we ignore the small devia- 



tions implied by Eq.(|28|), such that [A L>R (Z)) 00 = (g z LR ) SM ■ 1, [A LtR ) u = (gf) S M • 1 and 



[Al,r(Z)]22,33 = {9r)sm ■ I- To account for gauge boson mixing effects, we simply replace 



the first (zero order) terms in Eqs. (35) and (36) by [Al,r(Z)]q of Eqs. (29) and (30). 



We are now ready to compute the corrections to S(g1 e R )sM — from fermion KK mixing. 
Starting from Ai(Z)and using M i = and [Ai, tR (Z)]22,33 = (g^jsM • 1 we have: 



A L (Z) = [A L (Z)} 00 + {{g z R ) SM - (gf) SM )) 



M 02 M^ 2 Ml + M Q3 M^M\ : 



H)2 



03^33 



08 



(37) 



where (g R ) 



Ism — (gf)sM = 1/2. The block matrices M 02 and M 03 inherit their structure 
from the 5D Yukawa Lagrangian and are given by: 



Mi 



02 



X 



y e F e r e m x 



u 2 y T F e r^x 



(He + xlfy)F e r^x (w% + y e 2 f^)F e r^x (uy T + yifV}F e r^x 



(3* 



{y u + ei)Fir u im x e 2 Fir v wl x (e 3 + t x )F t r v im x \ 
Mm | e 3 FtrZoiX {y v + ei)Frt 01 x hF^x . (39) 

(e 2 + e x )F(r\ Q1 x e 3 F £ rl Ql x (y u + ei)F e r% 01 x ) 

We expect the most dominant contribution to come from the term proportional to 
Mq 3 M 33 M\ 3 induced by the interactions of the LH zero modes £® with the "custodians" 

Since t\ 2 3 come from the Z 3 preserving 5D operator k~'~ 7 ^ 2 £lH^ 2 i , r , 



s(i- + ) r,( 1 * + ) 



r 



1 R 



we expect them to be "rotated away" with V[ at 0{ei^,3, x -, fx) an d indeed using M ; 



-11 "f ■ 

1V1 KK 



83 



1 and Eq. (7), we factorize out {Fir\^ 2 v 2 j (M K ^ K ) ) 2 and obtain: 



(VlYM 03 M 3 - 2 M^Vl oc 



\y u \ 2 + lRe{y u e x ) - 2 fRe{y v e 



-^Re(yJ x ) 
-fRe(y u e x ) 



3 yyv^-x) 

y u \ 2 - fRe(&,e 



x) 



Re(y u e 



X) 



-^-Re(yJ x ) 
fRe(^e x ) 
\y v \ 2 - §Re(#„e x ) 



(40) 
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The above expression is already given in the mass basis and to get the contribution to 
the Zfxe coupling we simply take the (12) element of the matrix in Eq. (40), recalling the 
prefactors: 

,2 



KK — 



-UJ 



Re(y u e x )(F e r» l0l ) 2 



L KK 



9 X 10" 



(41) 



where we used r 



101 



1. We now proceed to obtain the 



0.72, F e ~ 0.1 and set \y u>ex \ 
correction coming from the M^M^ Mq 2 term in Eq. (37). We expect this term to generate 
a suppressed contribution due to the approximate alignment of M Q2 and M 00 . Notice that 
due to the degeneracy of q , we can write M 02 as 

M 02 = M 00 diag [F-^otTXr] , (42) 
using this compact form it is straightforward to transform this correction to the mass basis: 



(V^M 02 M^MlX = (V*)tMooV*(V*)tdiag 



diag(m eiM>T )(V^) t diag 



r, 



01 ) 



r e,fi,T' 00 1V1 KK 
„e,At,r\2 



/ 



i e,u,T\' 

(r r ) 



V 



r e,fj,,Ti 00 



/ 



n=l 



Vr diag(m eiAljT ) 



01 



y ^e,M,r r 00 



(43) 



/ 



It is now straight forward to estimate {5g^) KK by taking the (12) element from the above 
expression. The dominant term acquires the form: 



(Sg 



L IKK 



m e m M Af(ro\) 2 



F r M M m(1) 

/* 00 iw XX 



6 x 10 



-li 



(44) 



where we set all Yukawas to 1 in magnitude and use ~ 0.004, — 0.8 and Tq ~ 0.87. 
This concludes the calculation of KK fermion mixing contribution to Sg^ e . 

We now turn to the RH couplings, in which case we have only one term contributing, 
proportional to M\ M^ Mi and thus Eq. (36 simplify to: 



A R (Z) = [A R (Z)} 00 + ((g z L ) S M - (9 Z r)sm) 



(45) 



where (gf,)sM — (9r)sm = — 1/2, Mn = Mj£# • 1 and the various entries of M w are given by 
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(Ve + x[fl)F e r e 1Q x (Hp + xif^F^x {y T + x f 3 r)F T rl Q x 



M 



10 



y e F e r e w x 



u 2 y T F T rl x 



(46) 



% + x\f^)F e r e w x (u 2 y^ + y^f^F^x (uy T + ylfl)F T r T w x 



Once again, the degeneracy of q implies that M 10 is approximately aligned with M 00 and 
can be written in the following way: 



Al- 



io 



M 00 diag(ri M '7roo M,T )F; 



(47) 



The above equation simplifies the transformation of Eq. (45) to the mass basis and we get: 
(V^Ml M^M 10 V* = (Vfitdiag )v^M^V[ 



(VtfMnV&VtfdiBg 



10 



V, 



00 



R- 



(48) 



We use the same steps taken for the LH couplings in Eq. (43) to simplify the long 
expression in Eq. (48). The resulting expression for [AA R } mass = (T^MfoMxfMioVj reads, 



E (^)™(^WKJ 2 diag 



n,k,l=l 



10 

00 



(vS»(VSydiag 



10 
00 



(49) 

We realize that due to the degeneracy of q the LH diagonalization matrices are absent 
in the above expression, implying that the off diagonal elements of [Av4ijJ™ ass will come 
from terms proportional to the NLO corrections of Vr and will also consist of additional 
cancellation patterns, induced by the near degeneracy of the Higgs-flavon overlap correction 
factors tI'q qq. Specifying to the RH Zfie coupling, we extract the dominant terms and obtain 



KK 



m 2 A e 1 r 



10 



TO 



00 



10 
00 



_L0 

fi 



+ 



10 



-12 



(50) 



where we have omitted terms, which are at least 0(f x ~ 0.05) suppressed compared to those 
appearing in the above equation. To get the numerical estimation, we assigned all Yukawas 
to 1 and used r^o — 0.895, ~ 0.87 — 0.86, rf 
summarize our analytical estimations we have: 



0.754, r 



) ' 10 



0.74 and r[ — 



0.73. To 



9x10" 



(Sg 



R 



5x 10 



-li 



(5g£) KK « 10- 12 . (51) 



4.3 FCNC protection in the brane localized RS-A4 setup 

Since both the gauge and fermionic contributions, calculated in the previous sections, come 
from "cross-brane" effects it seems natural to believe that the brane localized RS-A 4 setup 
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is completely protected from all sources of tree level FCNC. This point was discussed in [H], 
yet the fermion KK mixing effects, were not taken into account there. Indeed, in the brane 
localized case we have = 1, which together with the degeneracy of q implies there can be 



no gauge boson mixing induced contributions to tree level FCNC (Eq. (30)). Similarly, the 
analytical approximations for {Sg1 e R ) KK (Eqs. (41), (44) and (49)) suggest that these sources 



of FCNC are turned off as well. Going to higher order in the expansion of Eqs. (35) and (36) 



won't do the trick since all blocks of the KK+zero mass matrix inherit their structure from 
the 5D Yukawa lagrangian, which is "form diagonal" in the brane localized case. Due to the 
Z% preserving VEV of $ (Eq. the interactions with the custodians (e^, r^), which 
mimic the Dirac interactions in the neutrino sector, are also "form diagonalizable" , namely, 
diagonalized by U(cu) (Eq. ([6])) and thus can not contribute to FCNC. 

The absence of tree level FCNC is not surprising since the LH zero modes are degenerate 
and the RH zero modes require no rotation to be transformed into the mass basis. This means 
that one can always find a basis in which the kinetic and mass terms are simultaneously 
diagonal, which simply forbids FCNC. Nevertheless, we have shown in the context of one 



loop contributions to dipole operators in the quark sector 10 , that the degeneracy between 



the KK masses, which is even stronger in the charged lepton sectors, implies that a non 
perturbative rotation of the degenerate KK blocks is needed before we can act with standard 
techniques of non-degenerate perturbation theory. Such a rotation, mixes the KK states 
maximally with each other (generation+BC) and might still induce small contributions to 
tree level FCNC. To verify the extent to which such a situation contributes to FCNC we shall 
attempt to analytically diagonalize the zero+first KK modes mass matrix. This seemingly 
hard task, shall be greatly simplified due to the magic properties of U(lj). We start by 
writing the 9x9 KK-KK block of the full mass matrix in terms of 3 x 3 building blocks, 



/ 
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9x9 
KK 



(X 
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KK 
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e,ti,T\ 







V 



ft 



KK) 



222 



diag(M XK ^ 




;(!+") 
KK 



■ 1 



(52) 



where v = 174 GeV is the Higgs VEV and Y„ is a dimensionless coupling matrix, inheriting 



its structure from the Dirac neutrino mass matrix (Eq. (13)) and is given by: 



Y 



KK 



Vv + ei e 2 e 3 
^3 Vv + ei e 2 
£2 £3 Vv + ei 



(53) 



Because of the symmetry between the e2, 3 entries, which is inherited from the Z% pre- 
serving VEV of <3>, we expect Yj CK to be diagonalized by U(u) and indeed we realize that: 
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Y 



KK 






ei + ue 2 + w 2 e 3 







y v + e x + u 2 e 2 + we 3 



■SyK K(diag . ) 



(54) 

Observing the above features of U(ui), we realize that the first trivial rotation we can 



impose on the KK mass matrix of Eq. (52), will be block diagonal, where each 3x3 building 



block is identified with either U (u or the identity matrix. Such a rotation will leave the Z 
coupling matrices, A LjR (Z) unchanged, since each 3x3 block in A L ^ R (Z) is proportional 
to the identity. In particular consider, the matrix Uj 2 ^ 12 = diag [Z7(a>)[l], 1, U(u), U(oS)\. 

Acting with jj^ 2x12 Q n M^j will make each 3x3 block generation diagonal, which means 
that the nine dimensional degenerate subspace, which had to be diagonalized has splitted 
into three, according to the "type" (BC) of KK modes. Assuming the degeneracy of KK 
masses, the rotated KK+zero mass matrix, m!^! is proportional to: 



/ diag(y, 
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diag(^ ejAtjT F ei 





t I 1 e,/i,T' 00 / 
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KK 



F t r v 222 



diag(y e ^ T )r e 1 f' T 
l/x 




)) 



kk t t r \m » 

\rv{diag.) p v 
KK 





l/x 



J 

(55) 

where x = v/M^k- We realize that all 3 x 3 blocks of the above matrix are diagonal, namely 
we can treat each blocks as a set of three complex numbers, without worrying about the 
inversion of parametric 3x3 matrices. Recalling the independence of U(u) from the LO 



Yukawa couplings, we realize that the 9x9 KK block of the matrix in Eq. (55) has the 
characteristic form: 



M 



Full 



KK 

(£/)) «M^e 



9x9 



l/x A B 
aA* l/x 
bB* l/x 



(56) 



where A = di&g(y e ^ T )rlf ,T , B 



Y£f»-) F tf u an d a(b) 



ii 



0.763, rft ~ 0.758, r T n ~ 0.753, r\ 



ii 



and r 



222 



0.19, we realize that a 



analytically by the transformation, M K £ 



ch.(diag.) 



diag^^^ / r^^j ) . Using 
0.749, r e 22 ~ 0.383 and r^ 2 ~ 0.367, r T 22 ~ 0.357 
1/2 • 1. The resulting can be diagonalized 



{V^ K )^Mf k V^ K , where are given by: 



KK 
L,R 



( 

\ B \ e ±i(e B -e A ) 

7w 2 



-|B| 2 



\A\ 



V2 

\A\e ±i( - e B- e A) 



2 + \B\ 2 ) 
\B\ 



V2 



2 +|SR) 

m 



(57) 



\ VlWP y/2(\A\* + \B\*) y/2(\A\* + \B\*) ) 
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and where M c ^ diag) = diag (l, 1 - 3/2a/|A| 2 + \B\ 2 , 1 + 3/2y/\A\ 2 + \B\ 2 ^J . We are now 

left only with the non vanishing zero-KK blocks- they are mixtures of the flavor diagonal 
blocks, diag(t/ ei/1)T F ei(UjT (£)rQ^Qs) and Y^ a9 '^ Fir\ m . Therefore, the remaining (standard) 

perturbative diagonalization of M^^ dia9 won't contribute to FCNC, to all orders. If we 
will consider the explicit values of the overlap correction factors, r nm above, flavor violating 
couplings may be generally generated at 0(v A /M^; K ), but will be further suppressed by 
quartic differences of the nearly degenerate r nm . Such contributions are negligibly small 
((9(1CT 18 )) and will thus not be further discussed. Finally, it is important to stress that the 
above diagonalization scheme can be easily generalized to include an arbitrary number of 
KK modes and will yield the same results in this case. 

4.4 The RS-A4 contributions to fi — > 3e and fi — > e conversion 

The effective Lagrangian for fi — > 3e deacy and \i — V e conversion in nuclei is written in terms 
of four fermion operators. These lepton flavor-violating Fermi operators are traditionally 
parametrized as p7] 

£ = -—^ [03 '^r){zr1^r) + 94{eL'fHL){e L 'fe L ) + g 5 (e R YfJ'R)(e L -f fl e L ) 

G 

+96(e L ^^ L )(e R -f ll e R )} + -j=e^{v - 075)// J^yy M (u 5 - a 9 7 5 )g, (58) 

v q 

where in the normalization we use, v q = T| — 2Q q sm 2 9. The axial coupling to quarks, a q , 
vanishes in the dominant contribution coming from coherent scattering off the nucleus. The 
#3,4,5,6 are responsible for fi — » 3e decay, while the v, a are responsible for \i — > e conversion 
in nuclei. The rates are given by (with the conversion rate normalized to the muon capture 
rate): 

Br(> -> 3e) = 2(g 2 3 + g\) + g\ + g\ , (59) 



p e E e G 2 F F 2 m la*Z A 
n 2 ZT 



Br(/, -+ e) = " ' :„.;,. " " e// Q> 2 + a 2 ), (60) 

capt 



where the parameters for the conversion depend on the nucleus and are calculated in the 



Feinberg- Weinberg approximation 38 and we write the charge for a nucleus with atomic 



number Z and neutron number iV as 

Q N = v u (2Z + N) + v d (2N + Z). (61) 
The most sensitive experimental constraint, Br(/i — > e) < 6 x 10~ 13 , comes from muon 



conversion in ||Ti 21 , for which 



10 6 

-E e ~p e ~m At , F p ~ 0.55, Z efr ~ 17.61, r capt ~ 2.6 . (62) 

s 
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Using these couplings one can estimate the coefficients of the 4- Fermi operators in ( 58 ) 



#3,4 



£5,6 = 25g£ : 6g% 



L,RJ 



(v ± a) 



25g 



(63) 



where Sg^ e R contains both the gauge and fermionic contributions to the Zfxe coupling. 



Recalling our collective estimations of the gauge and fermionic contributions (Eq. (51)) 



obtained for Yukawa couplings with magnitude 1 and random phase, we have (dg^^G+KK ~ 
8 x 1CT 8 (5 x 1CT 11 ). We expect the results of a numerical scan over the magnitudes and 
phases of the Yukawa couplings to be centered around these values if the effective field 
theory approach we used (Eq. (34)) is reliable. To find the exact result we perform a 



numerical diagonalization of M^ u « and use the resulting diagonalization matrices, ^ fulI ' ) 
to rotate the 12 x 12 coupling matrices Al,r(Z) = diag ([Ai ) K(Z)] o,n, 22,33) to the zero+KK 
mass basis, where we have (Sg^) mass = [A L:R (Z)}~ e(0) . Once we obtain (Sg^)^f a ~ s M we 
perform a matching to the operators contributing to /i — > e, 3e (Eq. ([58]) ) using Eqs. (59) 



and (60). Using the points satisfying the 3a bounds on the neutrino mixing angles from 
the scan of Sec. [3J we generate a sample of 40000 points, consisting of the e2,3,4, x , x e 2 3 
and y~2 3 parameters. The distribution of the values of each of these parameters, within the 
sample, is practically indistinguishable from the original distribution from which it was taken 
from. Namely, all parameters are still complex numbers with random phases and magnitudes 
normally distributed around 1 with standard deviation 0.5. Since BR(fi — > e, 3e) oc (of ^) 2 
(|60~1)) and 5gi generally dominates, we can plot the constraints coming from 



(Eqs. (J59J) 

the various cLFV experiments as vertical lines in the 5g^ z — $g R Z plane. The results are 
depicted in Fig. [7] (left). The results of a separate scan, independent of constraints from 
the neutrino sector, reveal negligible differences and are thus not discussed separately. By 
observing Eqs. pOl - ( 22 ) we recall that all mixing angles receive contributions from four 



NLO Yukawa couplings in the charged lepton sector £23 and y\ z . Consequently, there are 
simply too many input parameters of 0(1) and random phases, which can't be significantly 
constrained by the experimental data. By observing Fig. [7] we realize that the numerical 
results are in an excellent agreement with the estimations of Eq. (51). This strengthens 
our confidence in the effective field theory approach adopted from 
from SINDRUM II [21 



36 



The current bound 
around 35% of the cross-brane RS-A 4 setup and is 



"eliminates" 

easily and naturally satisfied for Air ~ 1.5 TeV. The correlation between the LH and RH 
Zfie couplings is attributed to the fact that they are both generated by cross brane effects, 
induced by v x (Eq. When cross brane effects are neglected, the Z/je couplings are 

negligible (right plot of Fig. [7]), while corrections to the TBM pattern can only come from 
higher order corrections to the heavy Majorana and Dirac mass matrices. Recall that in this 
case, achieving #13 ~ 6c/ \/2 is slightly less natural from the 5D theory point of view (requires 



y e4 \ ~ 6) (Sec. 3.1). We conclude that if no /i — > e, 3e events will be actually observed in 



Mu3e, MuSIC and Dee-Mee the "cross talk" RS-v4 4 model will be severely constrained and 
less appealing. On the other hand, the predictions for the brane localized realizations of 
RS-A 4 are extremely far from the reach of Mu2e (26], COMET [23], PRIME [24] and other 



future experiments. Most probably, the same situation will not hold true for the fj, — > ey 
decay induced at the one-loop level. The current upper bound is Br(/i — > cy) < 6 x 10 -13 [39] 
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Figure 7: The RS-A^ predictions for the anomalous LH and RH Zfie couplings in the presence 
(left) or absence (right) of cross brane interactions. Each point represents the contributions 
coming from both gauge boson mixing and KK fermion mixing, including all dominant higher 
order and cross talk effects. The dashed lines represent the maximum sensitivities of past, 
present and future LFV experiments taking place at FERMILAB and J-PARC. 



and is expected to improve by roughly two orders of magnitude in the next few years. This 
concludes our analysis of cLFV in RS-A 4 setups. 



5 Conclusions 

In this work we have studied the charged lepton and neutrino sectors of RS-A 4 (seesaw I) 
models excluding/including cross brane effects. For the neutrino sector, we have shown that 
cross-talk operators in the charged lepton sector and brane localized higher order correc- 
tions to the (heavy) Majorana and Dirac mass matrices induce significant deviations from 
TBM neutrino mixing, such that the experimental bounds can be rather easily satisfied. 
The mild differences between the normal and inverted hierarchy cases come from the terms 



proportional to £2,3,4 in Eqs. (20) -(22), associated with brane localized higher dimensional 
operators in the neutrino sector . The success rates for satisfying the 3<r bounds from all 
mixing angles simultaneously were shown to be £,nh{ih) — 10%(12.5%), which is slightly 
higher than those associated with "typical" A 4 seesaw models [31] and with anarchic mod- 



els 40 . Therefore, despite the fact that the recent measurements of 813 ~ Xc/V^ and the 



growing indications for the non maximality of #23 18 , 19 deviate significantly from the TBM 
pattern, models based on TBM mixing at LO with NLO corrections coming from both the 
neutrino and charged lepton sectors are still viable in explaining the neutrino mixing angles. 
The advantage of the RS-A 4 model remain in its relative simplicity and the fact that NP 
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contributions to cLFV processes are very suppressed. When specializing to the brane local- 
ized case, the resulting simplifications fix # 13 by a single (real) parameter and the remaining 
two parameters are first assumed to be equal up to a phase to look for possible correlations 
among the predictions. Indeed, we find a very strong correlation between the value of #13 
(close to 9°) and the the deviation from maximality of #23 towards the first (second) octant 
for the "positive" ( "negative" ) cases, as can be seen in Fig. [5j In the same context it is worth 
mentioning a slightly different approach 34 4l], in which the parameter space of the model is 
narrowed down to include only the VEVs of four flavons (3, 1, J/, 1"), corresponding to four 
input parameters in the Z 2 ^ preserving cases. These parameters are thus over-constrained 
by the existing neutrino measurements and the existence of a non trivial solution to this 
system of constraints reflects the suitability of A4 as a flavor symmetry for the lepton sector. 
Another interesting example is the so called "special" A4 models described in [31], where 
one is able to obtain a significantly better success rate £^ c ' ~ 50%. An example for such a 
model can be found in [42] . 

Turning back to the charged lepton sector, we have studied in detail the way in which 
KK mixing effects of gauge bosons and fermions generate flavor violating Z couplings in 
charge of tree level FCNC Due to the form diagonalizability of the LO mass matrices, the 
anomalous Z couplings are generated at the NLO by "cross-brane" effects which break A 4 
completely. We have obtained analytical approximations to the anomalous Zfie coupling 
using an effective field theory approach, in which all KK modes are integrated out. It was 
shown that the structure of the (non-custodian) corrections to the LH and RH Z coupling 
matrices (Eqs. (43) and (49)) are governed by the RH diagonalization matrices, implying a 
further suppression of f!^{jn e /m^) compared to the anarchic case and thus a partial protection 
against tree level FCNC. The results of the exact numerical diagonalization are in excellent 
agreement with the analytical approximations and are depicted in the left plot Fig. [7] We 
realize that the current upper bounds from SINDRUM II [21] eliminates roughly 35% of the 
model points for Ai R ~ 1.5 TeV, implied by EWPM. As for future Muon decay and conversion 



experiments, it can be seen that a non observation of // — > e conversion in DeeMEE 22 
(Br(/i — t- e) < (9(10 -14 ) already eliminates roughly 93% of the model points. To release this 
constraint we have to reduce the strength of cross brane interactions in the charged lepton 
sector, which forces us to re-match the neutrino mixing angles by gradually pushing the 5D 
Yukawas in the neutrino sector closer to their perturbativity limit (\Y\ < 10). A continued 



non observation of li -> e in MU2e [26], COMET [23], PRIME 24 and Project-X [27] will 
require very weak charged lepton "cross-brane" interactions and thus a too large value for 
Y. Turning off the "cross-talk" between the UV and IR brane, implies a complete protection 
from both gauge and fermionic KK mixing effects, as is the case for the earlier brane localized 
versions of RS-A4 [5j[6]. The latter protection stems from the fact that at LO no RH rotation 
is necessary (V£ = 1) and thus one can always find a basis in which the kinetic and mass terms 
are simultaneously diagonal, which completely forbids tree level FCNC. These conclusions 
also hold true once we consider the full KK mass matrix since all 3 x 3 blocks mimic the 
structure of the ZMA mass matrices. In particular, the blocks corresponding to interactions 
with the custodians (e^ R \ jj^ R , f^ R have the same structure as the Dirac mass matrix 
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(Eq. (13)), which in turn is also diagonalized by the same rotation Vf R = V[ = U(u). Thus, 
the RT>A.4 setup is interesting for cLFV, due to the non-custodial protection mechanism 
of the Z coupling matrices by virtue of the leading order "form diagonalizability" , which 
survives to a certain extent when cross-brane interactions are turned on. 

To summarize, we wish to stress that while anarchy is still viable in explaining the 
neutrino mixing pattern, we have chosen here the approach of looking for an underlying 
structure that can account for non trivial correlations among observables. It remains to be 
seen what future experimental data will tell us. Moreover, anarchic models generally imply a 
larger amount of flavor violation and are thus subject to stronger indirect constraints. Rather 
inevitably, the most common protection mechanisms developed to relax these constraints, 
need again to invoke a certain degree of underlying flavor symmetries. 
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